mu=1000;
f=@(t,x)[x(2);-mu*(x(1)^2-1)*x(2)-x(1)];
[t,x]=ode15s(f,[0,3000],[-1;1]);
plot(t,x)